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The  purpose  of  this  paper  is  to  present  a  new  local  explicit  method  for 
an  approximation  of  real-valued  functions  defined  on  intervals.  The  operators 
of  the  form  Qf  *  £  X  f  q  are  studied  under  a  uniform  mesh,  where  (q.  .  } 

i  1#K  1  |K 

comes  from  a  linear  combination  of  B-splines.  This  paper  contains  the 

definition  of  {q  },  comments  on  its  existence,  proof  of  reproduction  of 
i 

the  operator  Q  for  appropriate  classes  of  polynomials,  and  a  note  about  some 
applications. 
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SIGNIFICANCE  AND  EXPLANATION 


The  variation  diminishing  method  established  by  Schoenberg  and  the  quasi- 
interpolant  method  developed  by  de  Boor  and  Fix  take  the  form 

Qff  *  £  If  N,  where  (n  }  is  a  sequence  of  B-splines  and  {X  )  is  a 
1  1 / k  i/K  X 

sequence  of  linear  functionals.  This  form  is  convenient  in  practices.  We 
would  like  to  keep  this  form  but  replace  B-spline  ^  with  another  function 
qi  k'  we  consider  a  different  operator  Qf  ■  I  X^f  q^  where  qj^ 

has  small  support,  satisfies  q^^  ^(j)  *  ^ij'  and  *  f(x^).  Thus,  the 

operator  Q  becomes  interpolant,  and  Qf  is  in  a  class  of  the  so-called 
"many-knot"  splines.  The  paper  proves  that  Q  reproduces  appropriate  classes 
of  polynomials.  This  operator  can  be  used  to  fit  curves  or  surfaces. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  CLASS  OP  LOCAL  EXPLICIT  MANY-KNOT  SPLINE 
INTERPOLATION  SCHEMES 

D.  X.  Qi* 

As  is  well  known,  it  is  very  important  to  study  both  theory  and 
application  of  local  spline  approximation,  such  as  the  variation  diminishing 
method  established  by  Schoenberg,  the  quasi-interpolant  method  developed  by 
de  Boor  and  Fix  and  so  on.  Those  authors  studied  operators  of  the  form 
Qf  *  X  X.f  H.  ,  where  (n.  }  is  a  sequence  of  B-splines  and  {X  }  is 

^  1  1|K  if  K  1 

a  sequence  of  linear  functionals  (see  [1],  [2],  [3],  [4]). 

The  purpose  of  this  paper  is  to  present  a  new  method,  to  get  an 
approximation  of  real-valued  functions  defined  on  intervals.  In  this  method, 

I  use  (q  }  to  substitute  for  {n.  }  mentioned  above  as  a  basic 

ifX  if K 

function.  The  functions  q^  k  possess  the  following  characteristics:  (i) 
small  support  (it  makes  operators  of  the  form  Qf  “  X  X  f  q  local);  (ii) 

ifK 

qi,k(j)  “  Here  I  would  only  like  to  discuss  how  to  construct  the  basic 

functions  {q .  .  }  under  X  f  *  f ( x. ) . 

i  f  K  i  X 

Let  A  be  a  uniform  mesh:  a  *  Xq,  b  “  x^,  x^  -  Xq  +  ih  (i  *  0,1,..., N), 

and  additional  nodes  x_^,  x_2,...  and  x^.^,  xjj+2'***  *  Let  ^(A,k) 

denote  the  set  of  spline  functions  whose  knots  are  (x^,  x^  +  j).  Then 

8f  e  s  (A,k). 
p 

This  paper  contains  the  following  three  parts:  (i)  definition  of  a 

certain  basis  (q.  .  }  of  S  (A,k)  and  comments  on  its  existence,  (ii)  proof 
l,k  p 

that  Q  reproduces  approrpriate  classes  of  polynomials,  and  (iii)  a  note 
about  some  applications. 
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1.  Construction  of  {q  } 

1 


Let  be  Schoenberg's  centered  B-spline  of  order  k  on  a  uniform 

partition,  i.e., 

Mk<x)  ”  k[-  2  -x>+  , 

and  let  1  s*  {-  (k-2) , . . . ,k-2).  Then  the  functions 

M^i-*),  i  e  1 


are  B-splines  of  order  k  on  the  knot  sequence  Z  +  k/2,  hence  independent 
over  the  points  1/2  by  the  Schoenberg-Whitney  Theorem  (6]  since 
M^fi  -  i/2)  ^  0  for  i  e  I.  Consequently,  the  functions 


are  independent  over  I. 
Y  s°*  <Yi5iPl  80  that 


Mk( *- j/2 ) ,  j  e  i 

In  particular,  there  exists  exactly  one  choice  of 


I  YM  (.-j/2) 
jex  3 


(1.1) 


satisfies 


qk(1)  "  60i#  811  1  e  1  *  d.2) 

Note  that  *  Y^  by  uniqueness  and  symmetry  (which  can  be  used  to  simplify 

the  calculation  of  Y)  and  that 


1  -  I  qk(i)  -  l  l  Y.M.  (i  -  j/2) 

iei  iei  jei  3 


l  Y  (  l  M  (i  -  j/2) 

jei  3  iei 


l  Yj 

jei  J 


since  l  M.  (i  -  j/2)  *  £  K  (i  -  j/2)  -  1,  all  jei. 

iei  i 

Mow  we  define 


(1.3) 


qi,k(,)  qk<*”1)  * 

The  following  are  the  table  of  coefficients  Y  and  drawings  of  when 

k  -  2,3,4. 
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“2  -x  ■  u  l  i 


“3  -2  -1  10  1  2 


Figure  1 

D.  X.  Qi  (1975)  has  already  constructed  a  class  of  many-knot  spline 
interpolating  functions  for  solving  curve  fitting  problems  ([2],  [5]).  The 
main  difference  between  the  previous  study  and  the  present  one  is  in  their 
basic  function.  ^  that  appeared  in  [2J  and  [5]  is  not  the  same  as  q^. 

2.  The  interpolation  scheme  leaves  fixed 

In  this  section  I  want  to  prove  that  Q  reproduces  certain  polynomials 
I  will  use  the  symbols: 
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sym  (a  #a  s»  l  a  a  ...a  , 

V  2  *  <V*'V  1  2  W 

v  e  {l,2,...,k>,  *  v  (i  +  j)  , 


(0> 


8ym0(— )  >:  5JU#  «  1  , 


P(y)  ..  k-1  .  k-3  .  ^  k-1 .  .  k. 

:■  symy(i - — ,  i - — # ...,i  +  — ~)/(w) 


The  letters  P^  denote  the  set  or  linear  space  of  all  polynomials  of  order 
k,  i.e.,  of  degree  <  k. 

Lemma  (simple  consequence  of  Marsden's  identity  for  a  uniform  partition 

[4]) 

XW  «  l  5<»\(*-i).  X  e  [a,bj 

1  (2. 

P  *  0/1 f • • • f k* 1  • 


Theorem  1  Q|„  -  1. 

-  pk 

Proof  It  is  enough  to  prove 

xy  -  I  <i)Wq.  .<x),  x  e  Ia,b] 

j  1 

y  *  /k-i  . 


(2.2) 


Mow  we  use  induction  as  follows. 

Evidently  (2.2)  holds  for  y  ■  0.  Let  us  assume  (2.2)  holds  throughout 
y  -  0,1,..., ra-1.  We  will  prove  it  holds  for  y  ■  m. 

Notice  (1.1) 

q,  .  (x)  -  I  Y.M.  (x  +  i  -  i) 

i*x  jei  3  *  2 

and  by  lemma 

<x  +  ^)W  -  l  ^M\(x  +  \  -  i),  y  -  0,1,...,k-1  . 
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■**♦■'*■*■"  \+k& 


’•*  .  . 


& 

ii  > 


Therefore 


Pu(x)  t-  I  Y.(x  +  |)P  -  I  5!y)<li  v(x) 

W  jei  j  2  i  1  i>k 


Since  I  Y.  *  1# 


jei 


j 


P„(x)  -  I  Y.  (  I  (!!>xw"v(i)v) 


■  -i  V  U  '  V 

jei  J  v*o 


I  V,>*  l 

jei  3  \*i 


xy  +  l  <J)(  l  Y.(|)V)xW' 

v»i  jei  J 


l-V 


xW  +  I  (J>PV(0>X 
V*1 


By  induction  hypothesis  and  (2.3),  (2.4), 

fit  ,, 

.nu_  m-v 


U-v 


xm  -  P  ( x)  -  I  <v>V0)x 
V=1 


l  '1"’’ ,,„<«>  -  J,  <>v<°>  I  “’""V01’ 


I  a;m)  -  I  <>v<0)(i)~  )q4  k(x) 
i  v-1  ' 


Set 


m 


(ra)  f-  Am)  _  l  (”)Pv(0)im"V  . 

1  1  v-i 


Then,  from  (2.3)  and  ^( v)  * 

°j(°>  ■  l  5l3)’i,k(0)  '  Ei 


(J)  ,Y”i('  2 


k-1  k-1, 


However 


(2.3) 


(2.4) 


-5- 


The  last  identity  is  gotten  by  using  a  well  known  fact  about  elementary 
symmetric  function. 

From  Theorem  1,  we  can  get  a  result  about  approximation  order. 

Theorem  2  If  f  e  Ck+^[a,b],  then  R^  :*  f  -  Qf 
s) , 

■o,  “  max 

a+(k-1)h<x<b-(k-1)h 

s  *  0,1, ...,k  . 


|R^a)(x)|  -  0(hk+1'8) 


3.  Applications  in  CAGD 

By  convention,  let  denote  a  set  of  ordered  points  in  Rn.  We  hope 

to  get  a  curve  through  {p^K  It  is  known  that  people  in  Computer  Aided 
Geometric  (CAGD)  like  and  are  used  to  the  parametric  form.  So  the  curve,  as 
may  be  imagined,  can  be  represented  as  follows: 

Qklt)  =  |  qkU"j)Pj  *  <3.1) 

We  can  get  with  ease  from  this  representation  and  (1.1)  in  case  of  k  *  3,4: 


It  is  simple  and  useful  in  CAGD  that  the  interpolating  curve  is 


represented  by  a  matrix • 


(i)  Firstly,  we  consider  a  quadratic  many-knot  spline.  Let 


t  e  [0,  — ] .  We  can  find 


(q3<t+1),  g3(t),  q3(t-1),  q3< t-2)  =>  (t2,t,1)/j 


7  5 
4  4 


0  2 
1  0 


(3 


(t  ft# 1 )Mj  / 


and  with  the  help  of  symmetry 


1, 


(t2,t,i)M3(Pi_1,Pi,P1+1,P1+2)T,  t  e  [o,  •j] 


Q3(t) 


(<i-t)2,  1-t/i)M3(Pi+2,Pi+1,Pi,Pi_1)T/  te  1] 
(ii)  Secondly  we  consider  a  cubic  many-knot  spline.  Let  t  e 


[0, 


Then 


3  2  /  7 

<q4(t+2),q4(t+1),...,q4(t-3))  *  (t  ,t  ,t,1)  /  — 


11  II 

12  9 


_10 

9 


1  _1_ 
4  36 


2 

4 


3  5 

2  ”  2 


3 

2 


‘I  ° 


V 


12 


2 

3 


0 

1 


1 

12 


=*;  (t3,t2,t,1)M. 


and  with  the  help  of  symmetry 


Q4(t)  = 


(t3,t2,t,i)M4(p._2,Pi_1,...,Pi+3)T,  t  e  (o,  ±)  , 

((1-t)3,  (1-t)2,  1-t,  1)M4(P.+3,  Pi+2^  •  •  •  t  e  I2'  11  • 
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p. 


Figure  2 

As  the  parameter  t  increases  from  0  to  1,  the  segment  on  the  many- 
knot  interpolating  spline  curve  will  be  traversed  from  P^  to  (see 

Figure  2). 

If  we  want  to  get  raany-knot  spline  surfaces  when  the  points  {p  }  are 

i*  J 

given  (i  *  0,1,..., N;  j  ■  0,1,..., M),  we  could  represent  the  surface  as 
follows: 


Qk(u,w)  -  I  l  qlt<u-V>qk(jj-w)Pv^ 


V  U 

0  <  v  <  N,  0  <  w  <  M  , 

and  this  satisfies  Qk(i,j)  *  pif  j* 

The  representation  by  matrix  for  k  ■  3  is: 

(I)  Q3(u,w)  «  (u2,u,1)M3PM3(w2,w,1)T,  0  <  u,  w  <  ~  , 


P  P 

i-1/j-1  i-1,j 


i+2, j-1 


•  o  •  see 


i-1 , j+2 


i+2, j+2 


i+2 


(PV,U> 


j+2 


v-i-1,  U»j-1 


(II)  Q3<u,w)  «  ((1-u)2,1-u,1)M3PM3(w2,w,1)T,  |  <  u  <  1,  0  <  w  <  |  » 


p  ,  (p  l1”1  J+2 

'■rv,y;v»i+2,  y*  j-1  * 
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(Ill)  Q3(u,w)  =  (u  ,u,1)M3PM3((1-w)  ,  1-w,  1)  ,  0  <  u  <  -,  , 


P  "  ^Pv,y^\)=i-l ,  p=j+2  ' 


Q3(u,w)  *  (1-u)2,  1-u,  m3PM*(H-w)2.  1-w,  1)T,  ~  <  u  <  1,  •j<w<1  , 


P  =  (P  >i_1  j_1 

'  V, W ' v*i+2 ,  W-j+2 


Their  figures  are  shown  in  Figure  3. 


i+2 « j-1 


r  hi  /  iv 
i  /  ii  / 


i+2, j+2 


i-1, j-1 


i-1 , j+2 


Figure  3 


In  the  case  of  k  =  4  the  representation  and  figures  can  be  given  as 
follows: 

3  2  T  3  2  T  1 

(I)  Q4(u»w)  =  (u  ,u  ,u,1)M4PM4(w  ,w  ,w, 1 )  ,  0  <  u,  w  <  —  , 


P  -  (P  )1+3  j+3 

‘*V,U  V=i-2,  U= j-2  * 


(II) 


Q.(u,w)  -  ((1-u)  ,  (1-u)  ,  1-u,  1 )M.PM . (w 


T,  3  2 


P 


(P  )i_2 


j+3 

U+j-2 


I 


(III) 

Q4<u,w)  *  (u3,u2,u/1)M4PM^((1-w)3f(1-w)2,1-w,1)T,  0<u<~,  ~<w<1  , 


P 


(P  )1+3 
y 1  \*“i— 2 , 


j"2 

j+3  * 


(IV) 

Q4(u,vr)  *  ((1-u)3,(1-u)2,1-u,1)M4PM4((1-w)3,(1-w)2,1-w, 1)T,  ~  <  u,  w  <  1  , 


-  ^"2  *“2 
P  *  lPv,y’\>«i+3,  y-j+3  * 


Pi+3,j-2 


i+3,  j+3 
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